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ABSTRACT 

To understand the impact of low metallicities on giant molecular cloud (GMC) structure, we compare 
far infrared dust emission, CO emission, and dynamics in the star-forming complex N83 in the Wing of 
the Small Magellanic Cloud. Dust emission (measured by Spitzer as part of the S 3 MC and SAGE-SMC 
surveys) probes the total gas column independent of molecular line emission and traces shielding from 
photodissociating radiation. We calibrate a method to estimate the dust column using only the high- 
resolution Spitzer data and verify that dust traces the ISM in the H Fdominated region around N83. 
This allows us to resolve the relative structures of H2, dust, and CO within a giant molecular cloud 
complex, one of the first times such a measurement has been made in a low-metallicity galaxy. Our 
results support the hypothesis that CO is photodissociated while H 2 self-shields in the outer parts of 
low-metallicity GMCs, so that dust/self shielding is the primary factor determining the distribution 
of CO emission. Four pieces of evidence support this view. First, the CO-to-H2 conversion factor 
averaged over the whole cloud is very high 4-11 x 10 21 cm -2 (K km s -1 ) -1 , or 20-55 times the 
Galactic value. Second, the CO-to-H2 conversion factor varies across the complex, with its lowest 
(most nearly Galactic) values near the CO peaks. Third, bright CO emission is largely confined to 
regions of relatively high line-of-sight extinction, Ay > 2 mag, in agreement with PDR models and 
Galactic observations. Fourth, a simple model in which CO emerges from a smaller sphere nested 
inside a larger cloud can roughly relate the H2 masses measured from CO kinematics and dust. 
Subject headings: Galaxies: ISM — (galaxies:) Magellanic Clouds — infrared: galaxies — (ISM:) 
dust, extinction — ISM: clouds — stars: formation 



1. INTRODUCTION 

Most star formation takes place in giant molecular 
clouds (GMCs). A quantitative understanding of how lo- 
cal conditions affect the structure and evolution of these 
clouds is key to link conditions in the interstellar medium 
(ISM) to stellar output. Achieving such an understand- 
ing is unfortunately complicated by the fact that H2 does 
not readily emit under the conditions inside a typical 
GMC. Astronomers therefore rely on indirect tracers of 
H2, most commonly CO line emission and dust absorp- 
tion or emission. These tracers are also affected by envi- 
ronment, so that assessing the impact of local conditions 
on GMC structure requires disentangling the effect of 
these conditions on the adopted tracer from their effect 
on the underlying distribution of H2 . 

One way around this problem is to use several in- 
dependent methods to measure the structure of GMCs 
in extreme environments, inferring the state of H2 by 
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comparing the results. Here we apply this approach 
to an active star-forming region in the Small Magel- 
lanic Cloud (SMC). Using far infrared (FIR) emission 
measured by the S pitzer Survey of the SMC (S 3 MC 
iBolatto et all l200l and SAGE-SMC ("Surveying the 
Agents of a Galaxy's Evolution in the SMC" , Gordon et 
al. in prep.), we derive the distribution of dust in the re- 
gion. We compare this to maps of CO and H I lin e emis- 
sion (jBolatto et alj|2003t IStanimirovic et alJll999Q . Dust 
traces the total gas distribution — of which the atomic 
component is already known — and offers a probe of 
shielding from dissociating UV radiation. CO is the most 
common molecule after H2 (and the most commonly used 
tracer of molecular gas) ; understanding its relation to H2 
in extreme environments is a long-standing goal. The 
CO line also carries kinematic information that allows 
dynamical estimates of cloud masses. 

The SMC is of particular interest because the ISM in 
dwarf irregular galaxies like the SMC contrast sharply 
with that of the Milky W ay. They have low metallic- 
ities (e.g., iLee et all 120061) , correspondingly low dust - 
to-gas ratios fe.g.. Ilssa et al.|[l990|: IWalter et alJ [2007). 
and intense radiation fields (e.g., iMadden et al J l2006f > . 
These factors s hould affect the format i on and structure 
of GMCs (e.g . IMalonev fc Blackjll988t lElmegreenlll989t 
lMcKedH98alPapadopoulos et all l2002t ■ iPelupessy et al l 
2006). Unfortunately, it has proved extremely challeng- 
ing to unambiguously observe such effects because the 
inferred structure of GMCs depends sensitively on the 
method used to trace H 2 . 

Virial mass calculations reveal few differences between 
GMCs in dwarf galaxies and those in the Milky Way. 
In this approach, one uses molecular line emission to 
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measure the size and line width of a GMC. By as- 
suming a density profile and virial equilibrium, one 
can estimate the dynamical mass of the cloud indepen- 
dent of its luminosity. Recent studies find the ratio 
of virial mass to luminosity for GMCs in other galax- 
ies to be very similar to that observed in the Milky 
Way (IWalter et al.l l200lL l2002t i Rosqlowskv et all 12003: 



Bolatto et al J l2003t llsrael et all 120031 ; iLerov et allhoOB 
Blitz et all 120071 ; iBolatto et alj|20dl ). Further, the scal- 



ing relations among GMC s ize, line width, and luminos- 
ity found in the Milky W ay (|Larsonlll98lUSolomon et all 
Il987t iHever et al.ll2008h seem to approximately apply to 
res olved CO emission in other galaxies, even dwarf galax- 
ies (jBolatto et al.ll2008D . 

By contrast, observations of low metallicity galaxies 
that do not depend on molecular line emission consis- 
tently ^j[uggesj^lju-g£_j^3^^ CO 
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suggest large reservo irs or ri 2 untracea Dy 
Israellll997bt iMadden et all[l997t iPak et al]Tl99& 



i et all 12002: iGalliano et all 120031: IRubio etaO 



2004HLeroy et alJl2007t iBot et a l. 2007). The most com- 
mon manifestation of this is an "excess" at FIR or sub- 
millimeter wavelengths with the following sense: towards 
molecular peaks, there is more dust emission than one 
would expect given the gas column estimated from H I 
+ CO. llsrael (|1997bl ) treated the abundance of H2 as an 
unknown and used this excess to solve for the CO-to- 
H2 conversion factor. He found it to depend strongly on 
both metallicity and radiation field. 

These two sets of observations may be reconciled if 
CO is selectively photodis sociated in the outer parts of 
low- metallicity GMCs (e.g . iMalonev fc Blackll 19881 : llsraell 
1988; B olatto et all [19991 a scenario disc ussed specifi- 
cally for the S MC bv llsrael et all (|1986l ) and lRubio eFall 
(|199lLll993ah . This might be expected if H 2 readily self- 
shields while CO is shielded from photodissociating ra- 
diation mostly by dust, which is less abundant at low 
metallicities. In this case, CO emission would trace only 
the inner parts of low-metallicity GMCs. 

Observations of the Magellanic Clouds as part of the 
Swedish- ESO Submillimet er Telescope (SEST) Key Pro- 
gramme (jlsrael et al.lll993l ) support this idea : the surface 
brigh tness of CO is very low in the SMC (|Rubio et al.l 
119911 ): SMC clouds tend to be smaller than their Milky 
Way counterparts, w ith little associated diffuse emission 
(|Rubio et alj|1993at llsrael et al J 12003) ; and the ratio of 
13 CO to 12 CO emission is lower in the Magellanic Clouds 
than in the Galaxy, suggesting that clouds are more 
nearly optically thin ([Israel et all 12003) . 

The SEST results are mainly indirect evidence. What 
is still needed is a direct, resolved comparison between 
CO, dust, and H2. Because dust emission offers a tracer 
of the total gas dis tribution that is independent of molec- 
ular line emission ( Thronson "eTaT1[l987l H98l iThronsonl 
1 1 9 881: llsr aelll 1 99 7bh . it allows such a test. If GMCs at low 
metallicity include envelopes of CO-free H2, then the dis- 
tribution of dust (after subtracting the dust associated 
w ith H I) shou ld be ex tended relative to CO emission. 

ILerov et al.1 (|2007h attempted this measure- 
ment^ They combined S 3 MC with IRIS data 
(Mivillc-Desch ines &: Lagachd 120051 ) to derive the 
distribution of dust and comp a red th is to the NANTEN 
CO survey by iMizuno etai] (|2001h . They derived a 
distribution of H2 ~ 1.3 times more extended than that 
of CO, suggesting that half of the H2 in the SMC may lie 



in envelopes surrounding the CO peaks. The resolution 
of the CO and IRIS data limited this comparison to 
scales of > 45 pc. SMC GMCs ar e often much smaller 
than this (e.g. IRubio et al.l Il993al IMizuno et all 120011 ; 
llsrael et all 12003) . Therefore while this measurement 
indicated that SMC GMC complexes may be immersed 
in a sea of CO-free cold gas, it was not yet a true 
comparison of dust and CO on the scales of individual 
GMCs. 

Here, we focus on a single region, N83/N84 (hereafter 
simply N83). This isolated star-forming complex lies in 
the eastern Wing of the SMC an d harbors ~ 10% o f 
that galaxy's total CO luminosity (|Mizuno etaD'2001). 
Combining FIR, CO, and H I data we attempt to answer 
following questions: 

1. What is the CO-to-H2 conversion factor, Xqo (i- e -j 
the ratio of H2 column density to CO intensity 
along a line of sight) in this region? 

2. Is there evidence that CO is less abundant relative 
to H2 (i.e., that Xqo is higher or that there is H2 
without associated CO) in the outer parts of the 
cloud? 

3. Is the distribution of CO consistent with dust 
shielding playing a key role in its survival? 

4. Can dynamical masses measured from CO kine- 
matics be brought into agreement with H2 masses 
estimated from dust? What is the implied distri- 
bution of H2? 

To meet these goals, we first estimate the dust opti- 
cal depth at 160/im, ti 60 (||3|). We demonstrate that 
tiqo traces H I column density in the (assumed) H I- 
dominated ISM near N83, make a self-consistent deter- 
mination of the dust-to-gas ratio, and then combine T160 
with the measured H I column density to estimate the H2 
column density in the star forming region (Q. Finally, 
we combine the resulting maps of T160 and H2 with CO 
and H I data to answer the questions posed above ($5|). 

2. DATA 

We use FIR imaging from two Spitzer surveys. S 3 MC 
mapped 70 and 160/im emission from most active star 
forming regions in the SMC, including N 83. More re- 
cently, SAGE-SMC observed a much larger area, includ- 
ing the Magellanic Bridge and nearby emission-free re- 
gions. We use a combination of these data sets carried 
out by Gordon et al. (in prep.) that dramatically im- 
proves the quality of the 70/im image compared to S 3 MC 
alone, thus enabling this analysis. At 36" resolution, the 
noise (la) in the Spitzer maps is CT70 = 0.13 MJy ster -1 
(70/im) and ai6o = 0.6 MJy ster -1 (160/im) in the neigh- 
borhood of N83. 

We compare the Spitzer data to the IRIS 100/im image. 
IRIS is a re-processing of the IRAS data carried out by 
Miville-Desch enes fc Lagachd ([2005). These data have 
~ 4.3' resolution. 

iBolatto et all (|2003h used SEST to map CO J = 2 -» 1 
and J = 1 — > emission from N83. The half-power 
beam width of SEST was 23" (J = 2 — ► 1) and 45" 
(J = 1 — > 0). The maps that we use were convolved to 
lower resolution during reduction and have final angular 
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Fig. 1. — H I (left) and FIR, emission at 70 (middle) and 160 fim (right) in a two degree wid e fiel d centered on N 83. A thin black contour 
outlines the region where we can clearly distinguish FIR emission from the background (see 3231 ■ 



resolutions of 38" (J = 2 — > 1) and 55" (J = 1 — > 0). 
The noise in the velocity-integrated maps is somewhat 
position-dependent. Over regions with significant emis- 
sion lcr is typically 0.16 K km s -1 (CO J = 2 — ► 1) and 
22 K km s 1 (CO J = 1 ->■ 0). 

IStanimirovic et all (l999) imaged H I 21-cm line emis- 
sion across the whole SMC. These data have angular res- 
olution 98" and sensitivity sufficient to detect H I emis- 
sion along every line of sight within ~ 1° of N83. We 
cor rect for H I optical depth and self-absorption follow- 
ing [StanimirovicieF^n (|1999l . their Equation 6) ba sed on 
the H I absorption study by iDickev et all (|2000f ). The 
maximum correction factor near N83 is ~ 1.3. 

To subtract emission associated with the Milky Way 
from the FIR maps fi j2.ll), we use the Parkes map of 
Milky Way H I from lBruns et al.l (|2005h . Galactic H I is 
distinguished from SMC gas by its radial velocity. These 
data have a resolution of 14'. 

We move all data to three astrometric grids: one cov- 
ering the entire SMC, a two degree wide field surround- 
ing N 83 (Figured}, and the SEST field. In t he SEST 
field, we use the kernels of lGordon et al.l ()2008bl ) to place 
the 70/im image at the 160/im resolution (~ 36"), which 
matches that of the SEST CO J = 2 -> 1 data (38") 
well. We also convolve the 70 and 160/im maps to the 
55" resolution of the SEST CO J = 1 data. In the 
two degree field near N83, we use a Gaussian kernel to 
place the 70 and 160/im data at the 98" resolution of the 
H I. Over the whole SMC, we degrade the 70 and 160/im 
images to the 4.3' IRIS resolution. 

2.1. Additional Processing of the FIR Maps 

For consistency among the 70, 100, and 160 fira data, 
we move flux densities at 70 and 160 fxm from the MIPS 
scale (which assumes F v cx v 2 across the bandpass) to 
the IRAS scale (which assumes F v cx ^ -1 ). We do so 
by dividing the 70 map by 0.918 and the 160/im map by 
0.959. 

We subtract Milky Way foreground emission from the 
100 and 160/im maps. We estimate this from Galactic 
H I assuming the average c irrus dust properties measured 
by (|Boulanger et al.lll996[ K At 100/im we use their fit di- 
rectly; at 160/im we interpolate their fits assuming a typ- 
ical cirrus dust temperature (T = 17.5 K) and emissivity 
{(3 = 2). 



To refine the foreground subtraction, we assume that 
H I and infrared intensity from the SMC are correlated 
at a basic level. As the column density of SMC H I ap- 
proaches 0, we expect the IR intensity of the SMC to 
also approach 0. Therefore, we adjust the zero point of 
the IR maps using a fit of IR intensity to N (H I) SMC 
where iV(H I)smc < 2 x 10 21 cm -2 (we subtract the fit- 
ted y- intercept). This leads us to add 0.3 MJy ster -1 
at 70 fxm, subtract 4.4 MJy ster -1 at 160/im, and sub- 
tract 0.5 MJy ster -1 from the IRIS 100/im map. These 
offsets are a natural consequence of the uncertainty in 
the reduction and foreground subtraction (which must 
remove zodiacal light, Milky Way cirrus, and any cosmic 
infrared background). Deviations from the average cir- 
rus properties are particularl y common , being observed 
near a number of galaxies bv lBot et alJ (120091) . 

Based on carrying out this exercise in several different 
ways, we estimate the zero level of our maps to be un- 
certain by 0.25 MJy ster -1 at 70/im and 1 MJy ster -1 
at 160//m. We take these uncertainties into account in 
our calculations ( §3.2[) . To minimize their impact we 
only consider lines of sight with intensities well above 
the background, by which we mean /70 > 0.5 MJy ster -1 
and 7i6o > 2 MJy ster -1 after the foreground subtraction 
(i.e., twice the uncertainty in the background). 

2.2. A Word on Resolution 

In the rest of this paper we will combine the data de- 
scribed above in several ways. Two of these combinations 
lead to maps combining data with different resolutions. 
We comment on these here and the reader may wish to 
refer back to this section while reading the paper. 

First, we subtract a foreground component measured 
at 14' resolution from IR maps with 4.3' and ~ 36" 
(160/im) resolution. Any small scale variation in the 
Milky Way cirrus will therefore be left in our maps. This 
is only a concern in the diffuse region of the Wing (and 
so only in i j4.ll) . In N83 itself most lines of sight exhibit 
FIR intensities > 10 times higher than the foreground, 
so variations in the foreground are not a concern. 

Second, when estimating the distribution of H2 in N83, 
we derive the total amount of hydrogen (H I + H2) along 
a line of sight and then subtract the measured H I col- 
umn density. The total amount of hydrogen is based on 
FIR dust emission, measured at 36" resolution (or 55" 
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resolution when we compare to the SEST CO J = 1 — > 
map). The H I column density is measured at 98" reso- 
lution. We assume it to be smooth on smaller scales, an 
assumption born out to some degree by the reasonable 
correlation that we find between H2 and CO. Nonethe- 
less, the detailed distribution of H2 on scales less than 
98" (~ 29 pc) is somewhat uncertain. 

3. DUST TREATMENT 

We use the optical depth at 160/im, rieo, as a proxy 
for the amount of dust along a line of sight. For an 
optically thin population of grains with an equilibrium 
temperature Tdust, Ti60 is related to the measured 160/im 
intensity, I 160 , by 



Ti60 



'160 



B v (Tdust, 160/xm) 



(1) 



Here B v (Tdust > A) is the intensity of a blackbody of tem- 
perature Td us t at wavelength A. 

Calculating t\qq thus requires estimating Tdust- Be- 
cause only the 70 and 160/im maps have angular res- 
olution appropriate to compare with CO, we must do 
so using this combination. Unfortunately, T0/T6O does 
not trivially map to Tiust because the 70/im band in- 
cludes non-e q uilibr iu m emission f r om small grains (e.g., 
Desert et all 119901: iDraine fc ii 120071: iBernard et all 
20081) . We therefore take an indirect approach: we as- 
sume that most of the dust mass resides in large grains 
with equilibrium temperature Td us t that contribute all of 
the emission at 100/im and 160/im. We use T0/T6O to 
estimate Too / T.60 and then solve for Tj us t from 



Too 
T60 



100 
160 



-1.5 



By (T dust , 100/tm) 
B v (Tiust, 160/tm) 



(2) 



which assumes that dust has a wavelength-dependent 
emissivity such that t\ oc with j3 = 1.5. 

We derive the relationship between /70/T.6O and 
T00/T6O at the 4.3' resolution of IRIS, where both colors 
are known and exhibit a roughly 1-to-l relation. We then 
assume this relationship to apply to the smaller (~ 36") 
angular scales measured only by the Spitzer data. Near 
N83, the two colors are related by: 



^5° = o.24x 2 + 0.33a; + 0.45, where x = . (3) 

T.60 T60 

Note that this is not a general relation. It does not go 
through the origin and is only 1-to-l over a limited range 
of T0/T60; we n t an d apply over the range T0/T6O ~ 
0.15 - 1.2, where it is a good description of the SMC. 

3.1. Motivation 

In assuming that T00/T6O t races Td us t or its 
more sophisticate d analogs (e.g., iDale fe Heloul 120021 : 
IDraine fc lJl2007h . we follow several recent studies of the 
Magellanic Clouds dBot et alJ 120041: iLeroy et all 120071: 
IBernard et all 120081: [Gordon et alj|2008al ). ISchnee et all 
( 2005. 120061 . 120081 ) have demonstrated that a similar 
approach reproduces optical and near-IR extinction in 
Galactic molecular clouds, though with some systematic 
uncertainties. 




Fig. 2. — FIR color-color plot for the SMC. The x-axis shows 
To/^160, which is measured by Spitzer at high resolution but in- 
cludes contamination by non-equilibrium emission. The j/-axis 
shows /100/T6O, which is only available at 4.3' resolution but is 
more likely to trace exclusively equilibrium emission. Horizontal 
dashed lines show the temperatures associated with a few values of 
this color. Shaded contours show the distribution of data for the 
whole SMC; the lowest contour includes all data and the contour 
increment is a fact of 4 in data density. Black circles show mean 
and lcr scatter for data in a 2° field centered on N83 (binned by 
To/^l6o)- The histogram above the plot shows the distribution of 
over the SEST field (i.e., N83 itself). The dashed curve 
shows the color-color relation expected for a modified blackbody, 
which is not a good description of the SMC. On the other hand, the 
solid and dash-dotted lines, which show the relations that we use 
in our analysis, can reasonably predict -Z100/T6O from 

Figure [2] motivates our use of To / 7i60 (x-axis) to pre- 
dict T00/T6O (y-axis). Gray contours show the distri- 
bution of data for the whole SMC. Bins (filled circles) 
show data from a 2° square field centered on N83 (i.e., 
Figure [TJ. Both near N83 and over the whole SMC, the 
two colors show a reasonable correlation (rank correla- 
tion coefficient 0.7). 

Figure [21 also motivates our ad hoc treatment of the 
conversion between To / T60 and Too /T.60- A single mod- 
ified blackbody (the dashed line shows one with (3 = 1.5) 
cannot simultaneously describe the SMC at 70, 100, and 
160/im. The simplest explanation is that T00/T6O traces 
Tiust, while the 70/tm band includes substantial non- 
equilibrium e mission. We tested the possibility of using 
the models of IDraine fc Lil (|2007h . which include the ef- 
fects of stochastic heating, to directly derive dust masses 
from To/T60- However, the currently availab le "SMC" 
model s can not reproduce the data in Figure O iBot et all 
(|2004h and IBernard etaTI (120081) s howed that a simi- 
lar case holds for the iDesert et ail |l990) models. The 
main stumbling blo ck is repro ducing the observed 60/iin 
([Desert et all 1 19901 ) or 70/tm (jDraine fc Lil l2007f ) emis- 
sion. 

Equation [3] is not a unique description. A simple alter- 
native is a modified blackbody with twice the expected 
emission at 70/im. In this case: 



To 
T60 



2.0 x 



/ 70 
V 160 



-1.5 



By (Tiust, 70/im) 
By (Tiust, 160/im) 



(4) 



This is shown by the dash-dotted line in Figure [U It 
reproduces the data near N83 with about the same 
accuracy as Equation [3] If equilibrium emission sets 
Too / T.60 1 then Equation [5] implies that other processes 
(e.g., single-photon heating of small grains) contribute 
w 50% of the emission at 70/im near N83 (and across the 
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whole SMC). This is in reasonable agreement with the 
results for the Solar Neighborhood and several ne arby 
GMCs dDesert et al.lfl990t ISchnee et alJl^M l200l . 

The aim of this paper is not to investigate the details of 
small grain heating in the SMC, so we move forward us- 
ing our empirical fit (Equation[3]) . This appears as a solid 
line in Figured It is a good match to the data near N83, 
where the RMS scatter in the color of individual pixels 
about the fit is ps 0.04. In deriving uncertainties we use 
Equation[¥]as an equally valid alternative to Equation^ 

To convert from I\qq /Iieo to Td us t we assume that the 
SED along each line of sight is described by a modi- 
fied blackbody with t\ oc A - ^. At long wavelengths 
(A > 100/im), a blackbody spectrum with a wavelength- 
dependent emissivity is indeed a good descript i on of 
the integrated SEP of the SMC riAguirre et all 120031 : 
iWilke et all 12003: iLerov et alJl200l . We take (3 = 1.5, 
whic h is intermediate in the range of plausible values 
(e.g. , iDraine fc Led 119841 ) and a reasonable description 
of the integrated SMC SED from A - 100-1000/zm. This 
is not strongly preferred, and so we allow (3 from 1.0 to 
2.0 in our assessment of uncertainties. 

3.2. Uncertainties in tiqq 

We assess the uncertainty in t\§q by repeatedly adding 
realistic noise to our 70 and 160/im data and then deriv- 
ing T160 under varying assumptions. For each realization, 
we offset the observed 70 and 160/im maps by a random 
amount to reflect uncertainty in the background subtrac- 
tion; these offsets are drawn from normal distributions 
with la = 0.25 MJy ster -1 at 70jum and 1 MJy ster -1 
at 160/xm. We add normally distributed noise to each 
map. This noise has amplitude equal to the measured 
noise (SJ3]) and is correlated on scales of 36". 

We derive /100 /Iwo for each realization using either the 
polynomial fit (Equation[3]) or scaling the 70/im intensity 
(Equation d]), with equal probability of each. We add 
normally distributed noise to Iwo/Iwo with ler = 0.04 
(the RMS residual about Equations [3] and |4|) and then 
derive Td us t assuming (3 anywhere from 1.0 to 2.0 with 
equal probability. 

This entire process is repeated 1,000 times. We use 
the distribution of Monte Carlo T160S for each pixel to 
estimate a realistic uncertainty, finding individual mea- 
surements to be uncertain by « 40% (lc). We extend 
the same approach through our derivation of N (Hf IR ) 
in §4.41 In Appendix E] we discuss systematic effects 
that cannot be straightforwardly incorporated into this 
approach, two of which (blended dust populations and 
hidden cold dust) could impact tiqq. 

3.3. t 16 o and Extinction 

It will be useful to make an approximate assessment 
of the dust column in terms of F-band line-of-sight ex- 
tinction, Ay, and reddening, E(B — V). In the Solar 
Neighborhood E( B - V) = N (H) /5.8 x 10 21 cm" 2 
' Bohlin et al][l97l and r 160 = 2.44 x 10~ 25 cm 2 N (H I) 



A v [mag] = 1910 t 160 



(0) 



Boulanger et alj 119961 . studying the Galactic cirrus 
where we may safely assume that iV(H) « N(H I)). Then 

E(B - V) [mag] « 710 n 60 . (5) 

The reddening law in the SMC yield s Ry w 2.7 
(jBouchet et all 1 19851: [Gordon et al.ll2003f) . so that 



These equations assume the emissivity, t 1&0 /E{B — V), 
of Galactic H I but do not depend on the specific dust- 
to-gas ratio. 

Estimates of Ay and E(B — V) based on tiqq and 
Equations [5] and [6] agree well with opt ical- and UV- 
based measurements. I C apian et al.l ()1996[ ) compiled Ay 
for a number of SMC H II regions, including N83 and 
N84A/B (both of which lie within the SEST field). To- 
wards N83 they find Ay in the range 0.42-0.79 mag 
(mean 0.63 mag); towards N84A/B they found Ay from 
0.24-0.60 mag (mean 0.37 mag). Using their positions 
and aperture sizes, we derive Ay = 1.34 ± 0.36 mag 
and 0.93 ± 0.26 mag for the same regions. The optical 
and UV measurements are based on absorption toward 
sources inside the SMC. Therefore they will sample half 
the total line-of-sight extinction on average. Account- 
ing for this, our FIR-based extinction estimates are in 
excellent agreement with optical values. We find the 
same good ag reement for Sk 15 9, a B star near N83 
towar ds which iFitzpatrickl (|1984f ) and iTumlinson et al.l 
(2002]) measured E(B - V) w 0.05 mag, while we esti- 
mate E(B -V) = 0.08 ± 0.03 mag (see JO}. 

4. DUST AND GAS NEAR N83 

Following the method described in [J3j we calculate T160 
over every line of sight in a 2° field centered on N83 
(Figure [TJ and in the SEST field. In the process, we 
derive a median Td us t = 2 0.9 ± 1.5. This agrees with 
the T = 22 ± 2 K found bv lBot et all (|2004l ) for dust in 
the SMC Wing. The temperature in the N83 complex 
is somewhat higher, with median Td us t = 22.9 i 1.5 K 
and values up to ~ 28 ± 2 K. The hottest regions are 
coincident with the N83, N84A, and N84B H II regions. 

Our goal in this section is to combine T160 with the 
measured N(H I) to estimate N(H.2) via 

A(H™) = i(^-A(HI)) . (7) 
Here DGR is the dust-to-gas ratio defined by 



r 160 = DGR N(H) [cm" 2 ] 



(8) 



iV(H) = iV(Hl) + A(H 2 ), and Hf IR refers to the dis- 
tribution of H 2 derived using this approach. To calcu- 
late Hj 111 , we first compare tiqo and N (H I) in the area 
around N83 where the ISM is likely to be mostly H I 
( §4.1j) . This demonstrates that tiqq effectively traces the 
ISM and allows us to directly measure DGR in the dif- 
fuse ISM. We show that residuals about this tiqq-N (H I) 
relation come exclusively from regions of active star for- 
mation f H4.2p . We then adopt a reasonable value for the 
DGR in N83 itself and estimate N (H2) across the com- 
plex. 

4.1. H I and Dust Near N83 

In Figure [31 we plot T160 as a function of iV(H I) over 
the 2° field centered on N83. Most of the data are well- 
described by 
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Fig. 3. — Dust column, traced by rieo, (j/-axis) as a function of 
H I column density (x-axis) in a 2° field centered on N83. Black 
circles show average tiqq and la variation in bins 5 X 10 20 cm~ 2 
wide. The dashed line shows the median ratio Tieo/N(R I) = 1.4 X 
10 — 26 cm 2 . Dotted lines show the Hi column density for which 50 
and 75% of the pixels are well above the background (see §2.1; only 
such points are plotted). Contours show the distribution of data 
where point density is very high; the contour increment is a factor 
of two in data density. Dust and gas are reasonably related by a 
linear scaling over most of the field. The deviations to high Ti6o 
mostly coincide with N83 and other sites of active star formation. 

which is shown by the dashed line in FigureOU We expect 
that iV(H) ks N(H I) over most of this area. Thus, the 
clear, linear correlation in Figure [3] demonstrates that 
Ti6o traces the ISM well here and the slope is an estimate 
of the DGR in the diffuse ISM of the SMC Wing. 

Equation [5] i s consiste n t wit hin the uncertainties 
with results of iBot et ail (|2004h . who found ti 60 ~ 
(1.0 ±0.5) x 10" 26 cm 2 iV(Hl) [cm" 2 ] for the whole 
Wing (after adjusting for slight differences in Td us t, 
/3, and A). In the Solar Neighborhood, T160 ~ 
2.44 x 10 -25 cm 2 iV(H I) [cm~ 2 ] (jBoulanger et al.lll996fl . 
Comparing this to Equation[S]impUes that the DGR near 
N83 is 17^g° times smaller than the Galactic value. This 
agrees within the uncertainties with t he DGR found for 
the SMC Wing bv lLerovet all IpOOl . which is w I0±f 
lower than Galactic 10 . 

From Equations [5] and O we estimate N(R)/E(B — 
V) » 10±| x 10 22 cm" 2 mag" 1 . This matches the SMC- 
average N(R) /E(B - V) » 8.7 x 10 22 cm" 2 mag" 1 mea- 
sured by Fitzpatrick ( 1985) using IUE and confirmed by 
iTumlinsonet alj (j2002fl with FUSE. 

4.2. Residuals About the T160-H I Relation 

Equation [5] and Figure [3] demonstrate that a single 
DGR describes the region near N83 well. The notable 
exceptions are a small number of points with high T160 
relative to their H I column density. In Figure 2] we show 
the distribution of residuals about Equation [3J Contours 
indicate where our Monte Carlo uncertainty estimates 
yield 85, 98, and 99.9% confidence that the residuals are 
really greater than zero. 

10 Lcrov et al. (2007) made no correction for H I opacity. Doing 
so would improve the agreement with the present measurement. 



The neighboring panel shows the same confidence con- 
tours superimposed on an Ha image of the region near 
N83 (Winkler et al., private communication). The high- 
est residuals are associated with N83 itself. Other regions 
with higher-than-expected tiqq are also associated with 
concentrations of Ha emission. Ha emission indicates 
ongoing massive star formation, which in turn suggests 
the presence of H2. N83 also has significant C O emission, 
another signpost of H2 (|Mizuno et al.|[200lT ). If a large 
amount of the ISM is H2, we expect high residuals about 
Equation [9] even for a fixed DGR. 

4.3. The Dust-to-Gas Ratio in N83 

To derive from Equation [7] over the SEST field, 

we must know the DGR in N83 itself. We cannot mea- 
sure this directly because we do not have an indepen- 
dent measure of the H2 column. We might expect DGR 
in N83 to differ somewhat from that in the surrounding 
diffuse gas of the Wing: stars are more likely to form 
in regions with high DGR and the denser environment 
may shelter grains from destru ction by shocks or lead to 
gram gr owth (e.g.. lDwekj[l998D . In addition to our mea- 
surement of the diffuse ISM, we consider two pieces of 
evidence when adopting a DGR to use in N83: observa- 
tions of a nearby B star and the metallicity of the N84C 
H II region. 

FUSE and IUE Measurements of Sk 159: From FUSE 
and IUE absorption measurements, E(B — V), N(&2), 
and 7V(H I) are known towards Sk 159, a B0. 5 star near 
N83 (marked by a star in Figure |4]). H2 is detected 
but the column density is small (ss 2 x 10 19 cm~ 2 , 
I Andre et al.1 120041 ). The reddening as s ociated with the 
SMC is w 0.05 mag (|Fitzpatricklll984t iTumlinson et all 
2002) , though somewhat uncertain. The H I column 
measured from absor ption along the sam e line of sight 
is 2 ± 1 x 10 21 cm" 2 (jBouchet etal"1ll985| ). roughly half 
of the column inferred from 21 cm emission along the 
line of sight (two kinematically distinct H I components 
are visible in emission towards Sk 159; only one of them 
is seen in absorption, implying that Sk 159 sits between 
the two, behind the smaller one). These values imply 
N(B)/E(B—V) w 2-6 xlO 22 cm" 2 mag" 1 , or DGR w 2- 
7 x 10~ 26 cm 2 . 

Metallicity of N84C: iRussell fc Dopital (1990) mea- 
sured the nebular metallicity of the N84C H II region, 
which lies within the SEST field, finding 12 + log O/H = 
8.27, 2-3 times lower than the Solar Neighborhood value 
and among the highest for any region the SMC. Trans- 
lating metallicity into a DGR is not totally straightfor- 
ward, because the fraction of heavy elements tied up in 
dust may vary with environment. For a fixed fraction of 
heavy elements in dust, one would expect DGR oc Z . 
Fits to samples of galaxies yield power law relation- 
ships (DGR L j2£_Z2)_with = 1-2 
(e.g JLisenfeld fc Ferraralll998l : lDraine et al.ll2007T ). This 
would imply N(R)/E(B -V) ~ 2-7 x 10 22 cm" 2 mag" 1 
or DGR ~ 2-7 x 10~ 26 cm 2 . 

H I and rigo : Equation [9] offers a lower bound on 
the DGR — N83 is extremely unlikely to have a lower 
DGR than the surrounding medium (N(R)/E(B -V) ps 
10 x 10 22 cm -2 mag -1 ) and from absorption work we 
know that there is not a pervasive massive molecular 
component in the SMC. The magnitude of the resid- 
uals about this equation towards N83 itself also offer 
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Fig. 4. — (left) Residuals about Equation [9] the average relationship between rieo and H I in a 2° field centered on N83. The contours 
show where our Monte Carlo uncertainty estimate yields 85, 98, and 99.9% confidence that the residual is above zero. A thin gray line 
shows where we clearly distinguish FIR emission from the background (see S|2,U and the star indicates Sk 159, a B star observed by FUSE 
and IUE. Most of the region is well— described by Equation [9] but N83 itself shows higher Tigo than expected from H I and Equation [9] 
(right) The same contours plotted on top of Ho emission near N83. Regions with high rigo residuals are associated with sites of recent 
high-mass star formation. 



a weak upper bound on the quantity. If we assume 
DGR much above 3 times the value in Equation [5] 
then some lines of sight inside the SEST field will have 
significantly negative residuals. If the star- forming re- 
gion itself is described by a single DGR, then it must 
be roughly bounded by this value, which translates to 
N(R)/E(B - V) ~ 3 x 10 22 cm' 2 mag" 1 . 

Assumed DGR in N83: The relatively high metallic- 
ity and the measurement towards Sk 159 are balanced 
against our observations of a very low DGR in the 
nearby ISM and the requirement that Sh2 R n °t signif- 
icantly and systematically negative. The former suggest 
N{R)/E{B - V) ~ 2-7 x 10 22 cm" 2 mag" 1 , while the 
latter yields N(R)/E(B-V) ~ 3-10xl0 22 cur 2 mag" 1 . 
In the remainder of this paper we adopt assume that in 
N83 itself N(R)/E(B-V) ~ 5xl0 22 cm" 2 mag" 1 , which 
is intermediate in this range. Then 



Ti 60 = 2.8 x 1CT 26 cm 2 N(B) [cm -2 ] 



(10) 



This is twice the value found in the diffuse gas of the SMC 
Wing (Equation [9]) and more similar t o that found in the 
actively sta r-forming SMC Bar (e.g.. IWilke et al.ll2004 
iLerov et al.ll2007lh It is roughly consistent with observa- 
tions of Sk 159 and the metallicity of N84C. This DGR 
also leads to reasonable agreement between dynamical 
and dust masses in the star-forming region (i i5.4p . which 
was a factor in settling on this value. In Appendix |A1 we 
illustrate the effects of changing this value on our analy- 



4.4. H FIR in N83 

Combining Equations 171 and flOl we estimate ./V(H 2 IR ) 
from Ti6o and N(R I). From A(H FIR ), we calculate the 
molecular gas surface density, 



v fir 

Zj H2 



[M 



PC 



A(H™) 



4.6 x 10 19 [cm" 2 ] 



(11) 



which includes a factor of 1.36 to account for helium 11 
(after I Wilson et al.lll988l) . At the same time we estimate 
the extinction along each line of sight using Equation [6l 
Carrying out these calculations, we work with N(H I) 
only in average, because the resolution of the 160 /Mil 
and CO J = 2 -> 1 data are ~ 38", while that of the H I 
map is 98" f ^]) . 

The right column in Figure [5] shows the resulting maps 
of in N83 at the resolution of the SEST CO 2 -> 1 
(top) and 1 — > (bottom) data. The left column shows 
the CO maps. Note that the stretch on the H FIR images 
runs linearly from Eg 1 ^ = 100 M Q pc -2 to 500 M pc" 2 . 

Several systematic uncertainties may affect N (H 2 IR ) 
but are hard to quantify and so not reflected in our Monte 
Carlo estimate of the uncertainties. We consider these 
in Appendix fA) finding no strong reason to doubt that 
Equation [7] yields an approximate estimate of ./V(H 2 ). 



5. H^ 111 , CO, DUST, AND DYNAMICS 
5.1. Hf IR and H I 

Before we consider the relationship between CO, H FIR , 
and dust within N83, we briefly examine the transition 
fr om atomic (H I) to mo lecular (H 2 ) gas in the complex. 

iKrumholz et al.l (|2009f ) recently considered the transi- 
tion from H I to H 2 in galaxies. They argue that inside 
a complex of mixed atomic and molecular gas, the ratio 
of H 2 to H I along a line of sight (i?H2 = Sh2/^hi) is 
mainly a function of two factors: total gas surface den- 
sity (Shi + Sh2) and metallicity. Their calculations agree 
well with a variety of observations, including FIR-based 
estimates of Sh2 in the SMC at lower resolution. 

Comparing H I and H 2 in the area around N83, we in- 
deed observe a clear relationship between i?H2 and the 

11 In the rest of the paper, £^2 R includes this correction for 
helium, while A^(H2) or A^HJp^) refer to column density of H2 
alone 
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Fig. 5. — CO emission (left panels) and S]^ 11 estimated using Equation [7] (right panels). The top panels show CO J - 



and H^ IH at 



38" resolution (but see [[272] regarding the resolution of the H2 IR map). The bottom panels show CO J = 1 



1 emission 
emission 



and Hj 111 at 55" resolution. Dotted contours show the boundaries of the SEST map. In the CO maps, contours show Iqq from 1 to 
8 K km s — 1 spaced by 1 K km s — 1 . In the HJp 11 maps, contours indicate S^2 R from 100 to 500 Mq pc - 2 spaced by 50 Mq pc - 2 . 



total gas surface density. We show this in Figure plot- 
ting i?H2 against Shi + ^H2 R over the whole area where 
> °- We work at the 98 " ( 29 P c ) resolution of the 
H I map, with each point in the plot showing an indepen- 
dent measurement. For this analysis, we are interested 
in the gas associated with the star-forming complex it- 
self (not unassociated gas in front of and behind it along 
the line of sight). To remove H I unassociated with N83 
itself from Shi, we subtract the median Shi measured 
over the area shown in Figure Q] (53 M Q pc -2 ) from the 



measured Shi before plotting. This is only an issue for 
H I; Hf IR does not extend beyond the N83 complex. 

We overplot th e relationship between i?H2 and Shi + 
Sh2 predicted bv lKrumholz et all (|2009r ) for three metal- 
licities: Z = 0.5, 0.33, and 0.125 ti mes solar. Our 
data a re consistent with the shape of the lKrumholz etahl 
(2009) calculation. We find i? H2 = 1 at S H2 + S H i = 
68 ± 12 Mq pc -2 , which agrees well with their calcula- 
tions for Z 2-3 times lower than the solar value. This 
is roughly the metallicity measured for the N84C H II 
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Fig. 6. — Ratio of molecular to atomic gas, Ru2 = >l[j!> /^'lll : 
as a function of total gas column density Shi + Sh2 in the N83 
complex (after removing H I not associated with the complex). 
Each point shows an independent measurement at 29 pc resolution. 
The dotted line shows the sensitivity of our Hj IR map. We plot the 
th eoretical relation s hips b etween i?H2 and Shi + Sh2 calculated 
by IKrumholz et al. (2009) for several metallicities. 



region (jRussell fc Dopital fl990). However, it is signifi- 
cantly higher than the DGR that we ad opt (^4.31). which 
is clos er to the lower value. Because IKrumholz et al.l 
(2009) assume a linear scaling between dust opacity and 
metallicity when deriving these curves, this means that 
there remains some disagreement between our measure- 
ments and their results. Nonetheless, there is good qual- 
itative agreement in the shape of the curve and the fact 
that in N83 i?H2 = 1 at a significantly higher value of 
Shi + £r~2 t nan in a solar metallicity cloud. 

5.2. CO and Hf IR 

Figure [5] shows that the distributions of H FIR and 
CO share the same peaks and basic morphology. How- 
ever, the values of ico i n N83 are low compared to a 
Galactic molecular cloud, which usually show Ico ~ 
10 K km s 1 over a larg e area, not merely the peaks 
(e.g., iWilson et al.l 120051 ). By contrast, the values of 
(mean 180 M Q pc -2 ) are similar to the surface 
density of an average Galactic GMC ~ 1 20-170 M Q pc~ 2 
(jSolomon et alj|1987t IBever et al.ll2008fh 

This means that CO is faint compared to H F IR in N83. 
Over the SEST field X co is 



(X& 1 ) = 6.7±HxW 
<Xfe°>=7.9±^xl0 21 



cm 



K km s 1 

cm 
K km s _1 



(12) 



These ratios are 34 
tic conversion factor, 



and 40 times the Galac- 
taken to be Xg a i ~ 2 x 
10 20 cm~ 2 (K km B- 1 )- 1 (e.g., iStrong fc Mattel [1996L 
iDame et al.l l200lT ). This value agrees reasonably with 
previous FIR-based determinations of Xco in the SMC 
and N83: comp aring IRAS an d CO at selected point- 
ings in the SMC lIsraell (|1997bD derived X co ~ 67 X Ga i. 
Applying the same methodology to N83, iBolatto et al.1 



(2003) found X CO ~ 100 ± 50 X GaX . iLerov et all (p007h 
derived X CO ~ 50 A Gal comparing NANTEN CO, IRIS 
100/im and Spitzer 160 /im towards N83 (removing their 
correction for extent). 

The left panel in Figure [7] compares H FIR and Ico for 
individual lines of sight. We plot S™ as a function 
of Ico over the SEST field. We regrid the data so that 
each point corresponds to an approximately independent 
measurement over a ~ 10 pc (CO J = 2 — > 1) or ~ 17 pc 
(CO J = 1 — ► 0) wide box. Gray curves show fixed CO- 
to-H2 conversion factors, starting with Galactic (lowest) 
and increasing by factors of 3.33. 

As with Figure Figure [7] shows that despite the very 
low ratio of CO to H FIR , the two exhibit an overall cor- 
respondence. High Ico coincides with high S H 2 R and 
the reverse, so that a rank correlation coefficient of 0.7 
relates the two over the SEST field. 

The relationship between Ico and S^ : 2 R does not go 
through the origin. Instead, Ico = corresponds to 
roughly l£ IR = 50-150 M Q pc" 2 . This suggests the 
presence of an envelope of Hf IR with very little or no as- 
sociated CO. Unfortunately, this result is very sensitive 
to the adopted DGR (£ 14.31 and Appendix |A|) . If we take 
DGR at the upper end of the plausible range, the data 
are consistent with no CO-free envelope although CO 
emission is still faint relative to Sh2 R i n ^ ne SEST field. 
If we take DGR at the value derived in the nearby diffuse 
ISM, the surface density of the envelope is even higher 
~ 200-400 Mq pc -2 . Although the observation towards 
Sk 159 does not actually intersect the envelope in the 
latter case, it is very nearby and the low iV(H2) derived 
from absorption towards this star offers some circumstan- 
tial evidence against a very massive extended envelope. 

The other notable feature of this plot is that at very 
high CO intensity increases dramatically (the turn 
to the right at the top of the plot). We see this in both 
CO transitions, but the effect is more pronounced at the 
higher resolution of the CO J = 2 — > 1 data, suggesting 
that the bright CO-emitting structures are still relatively 
small compared to the SEST beam. The result is that the 
line-of-sight integrated ratio of H FIR to CO is lower for 
the regions of brightest CO emission, dropping to ~ 15 
times the Galactic value. Care must be taken interpret- 
ing these ratios because H FIR and CO emission almost 
certainly trace different volumes ( £|5.4[) . 

5.3. CO and Extinction 

In £JTJ we highlighted the role of dust in shielding 
CO from dissociating radiation. This may provide a 
simple expla nation for the upturn in CO intensity at 
high Sg IR . iLeaueux et all (|l99"l modeled CO emis- 
sion in SMC molecular clouds. For their typical cloud 
(njj ~ 10 4 cm -3 , illuminated by a radiation field 10 
times the local interstellar radiation field), they found 
that most CO emission comes from a relatively narrow 
region of the cloud centered on Ay ~ 1 mag. Outside 
this regime CO intensity is very weak, a scenario that 
qualitatively mat ches what we se e in the left panel of 
Figure E| (see also lBell et aTll2006D . 

In the right panel of Figure [7] we plot CO intensity 
as a function of line-of-sight extinction, Ay. We esti- 
mate Ay from Ti6o using Equation [S] For comparison, 
we mark Ay ~ 2 mag, the line-of-sight extinction that 
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Fig. 7. — The structure of H2, CO, and dust in N83. (left) Hip^ surface density, S^ 1 ^, as a function CO intensity, Ico (x-axis). Dotted 
gray lines show CO-to-H2 conversion factors of 1, 3.33, 10 ... 333 times the Galactic value, (right) Iqq (j/-axis) as a function of line-of-sight 
extinction, Ay (i-axis), estimated f rom Ti60 . The vertica l line shows the line-of-sight extinction from which most CO emission emerges in 
models of SMC molecular clouds by Lequ eux et alj 11994ft . The gray curve shows the relationship between Ico an d extinction observed in 
the Pipe Nebula (Milky Way) bv lLombardi et alj 112006ft , In both plots, each data point represents an independent line of sight. We show 
results for the J = 2 — > 1 transition (38" resolution) in black and the J = 1 — > transition (55" resolution) in gray. 

tion and it certainly seems likely that the radiation field 
incident on the H2 in N83 is much more intense than in 
the relatively quiescent Pipe. 

Small differences should not overshadow the similar- 
ities between the CO-extinction relation in the Milky 
Way and that in the SMC. Compared to the left panel in 
Figure [3 the right panel actually shows a striking sim- 
ilarity between Galactic and SMC clouds. We derive a 
CO-H2 conversion that differs with the Milky Way by 
a factor of ~ 30, while the relationship between extinc- 
tion and CO is only slightly offset. Figure [7] supports 
the hypothesis that shielding, rather than the distribu- 
tion of H2, determines the location of bright CO emis- 
sion. Here "shielding" refers to a combination of dust 
and self-shielding. Both processes are important to set- 
ting the location a t which most C is tied up in CO (e.g., 
Wolfi reTet al.lll993h and the effective shielding from both 
sources will be weaker in the SMC than in the Galaxy 
due to the decreased metallicity. 

Extinctio n may a l so be critical to a cloud's ability to 
form stars. iMcKed (|1989| ) proposed that ionization by 
an external radiation field plays an important role in set- 
ting cloud structure because it determines the degree of 
magnetic support. He predicted that clouds forming low- 
mass stars in equilibrium will self-regulate to achieve in- 
tegrated line-of-sight extinctions Ay ~ 4-8 mag. These 
extinctions are higher than the Ay ~ 2-3 mag that we 
find towards the CO peaks N83 or the average extinc- 
tion over the region, Ay ~ 1.5 mag. We can safely con- 
clude that the N83 region as a whole does not resemble 
th e equilibrium low-mass star forming cloud described 
bv iMcK cc (1989). If these equilibrium structures do ex- 
ist in this re gion, they must be compact relative to our 
10 pc beam. iBolatto et al.1 ([2008) find that the dynamics 
of CO emission in the SMC also appear to disagree with 
the predictions of [McKei ( jl989i ) but present several im- 
portant caveats to the comparison. The most important 



rough ly matches the depth from which iLeaueux et al.l 
(1994) predict most CO emission to emerge (see their 
Figures 2 and 6). They model a slab illuminated from 
one side while we estimate the total extinction along the 
line of sight through the cloud. Therefore Ay — 1 mag 
for them corresponds to Ay ~ 2 mag for us (though 
the actual geometry is likely to be much more compli- 
cated). We also plot the relationship between extinc- 
tion and CO intensity meas ured in the Pipe Nebul a (a 
nearby Milky Way cloud) by lLombardi et al.l (|2006t see 
their Figure 22). We convert Ak into Ay using their 
adopted Ay — Ak/0A12. They measure a scatter of 
roughly 2 K km s _1 about this relation. 

In agreement with ILeaueux et "all (|1994f ). we find that 
lines of sight with bright CO emission occur almost ex- 
clusively above Ay ~ 2 mag. Our maps lack the dy- 
namic range in Ay to test whether Iqo is indeed more or 
less independent of e xtinction well above this threshold 
(as m the Milky Way. lLombardi et al.ll2006tlPineda et all 
2008). In fact, Figure 3a of ILeaueux et al.l (|1994l ) seems 
a close match to what we observe: a shallow slope that 
steepens sharply around Ay of 2 mag (for us). The radi- 
ation field that they assume, 10 times the Galactic value 
is a rough match to what one would infer comparing Tdust 
in N83 (median ~ 23 K, max ~ 28 K) to that of Galactic 
cirrus (17.5 K) — median ~ 5, maximum ~ 15 12 — es- 
pecially when one recalls that this is integrated over the 
whole line of sight rather than tracing the radiation field 
incident on the cloud surface. 

N83 shows somewhat less CO at a given extinction 
than the Pipe Nebula. This is also in agreement with 
the models by lLequeux et al.l (|1994l ). which predict that 
CO from Milky Way clouds emerges from a broader range 
of Ay and lower values of Ay than in the SMC. They 
attribute the difference to lower rates of photodissocia- 

12 For our adopted /3 = 1.5, the magnitude of the radiation field 
heating the dust is roughly oc T 5,5 . 
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of these here is that iMcKed ()1989t ) explicitly consider 
clouds forming only low-mass stars, while N83 is quite 
obviously actively producing high mass stars. 

We emphasize that this comparison between CO and 
Ay is fairly robust. It does not depend on our choice of 
DGR, only on the adopted FIR emissivity (tfir/Av) and 
reddening law. The most likely biases in the emissivity 
(e.g., coagulation of small grains) will lower Ay, bringing 
our results into even closer agreement with those in the 
Milky Way. 

5.4. Hf IR and Dynamical Mass Estimates 

CO line emission also offers kinematic information. 
This is the basis of the virial mass method commonly 
used to estimate the masses of molecular clouds and 
derive CO-to-H2 conversion factors (e.g . IRubio et al.l 
H993 afc iMjjo^liaaa lArimoto et al1ll996T ) including in 
N83 (Bola tto et aH 2003: Isr ael et alf 2003: B olatto et all 
2008) . The potential pitfall of this approach may be seen 
from §5.31 if CO emission is confined to regions with ex- 
tinction above a certain threshold and these regions rep- 
resent only a fraction of the whole cloud, then velocity 
dispersion and size measured from CO observations will 
only be lower than their true values. Mass outside the re- 
gion of CO emission may exert pressure on the surface of 
the "CO cloud," but it is not straightforward to estimate 
the total mass of a cloud from observing kinematics from 
only part of it. As a result, in a low-metallicity cloud like 
N83, we expect virial masses from CO observations to be 
smaller than Hf IR (even over matched areas) because the 
latter also traces the outer (CO-free) part of the cloud, 
which exists in front of and behind the CO-emitting re- 
gion even over matched lines of sight. 

In N83, we have the advantage of an independent mea- 
surement of H FIR and observing the CO emission over 
a range of scales. Here we test whether these observa- 
tions can be reconciled using a simple model in which 
CO emission comes from only the inner part of a larger 
H2 cloud (as appears to be the case in N83). We con- 
sider a spherical cloud with a radially declining density, 
such that p oc r~ a , and a radius R beyond which p = 0, 
i.e., the model usual ly adopted (with a = 1) to calculate 
cloud virial masses (Solo mon etallfl98l 13 . We assume 
that the dynamical mass estimated from CO line data 
traces the mass of a fraction of this cloud, out to ra- 
dius rf rac . The ratio of dynamical mass to Hip 11 over a 
matched area, M^/M^ 1 is then a function of a and 
the ratio of the true radius of the cloud to the radius of 
the area being considered, r{ rac /R. The top left panel of 
Figure [H] shows this ratio for models with a from to 
2.0. 

To compare our observations to this model, we mea- 
sure the line width and radius of CO emission over a 
series of scales in N83. We consider intensity contours 
in position-position-vclocity space, beginning with the 
bright northwestern region and including progressively 
more of the cloud (but always including that region, see 
Figure [5]). We estimate the radius and line width of each 
region from the area (for the radius) and second moment 
(for the line width). To account for the finite resolution 
of SEST, the radius of each cloud is adjusted by 

13 We cap the density at its maximum value over the inner 3% 
of the cloud to avoid divergence. 



* = " ^ ■ (13) 

Here A c i ou( j is the area of the cloud and -Bbeam = 
OMF WHM is the "radius" of the beam ([Solomon et all 
1987;. We combine the RMS line width, er„, and cloud 
radius, R, to derive the virial mass via 

Af vir = 1040i^ [M ] , (14) 

with <j v in km s _1 and R in pc. For details of measuring 
the properties of extra galactic GMCs from CO emi ssion, 
we refer the reader to IRosolowskv fc Lerovl (|2006l ) and 
references therein. 

For each contour, we measure M v - lr / 'M F 2 R . We com- 
pare this ratio as a function of R to a range of density 
profiles and cloud radii. The resulting distribution of 
reduced \ 2 lS shown in the bottom left panel of Figure 
[8l Our measurements, along with the best-fit model are 
shown in the bottom right panel of the same figure. 

The best-fit model has p oc r~°- 6 and R = 70 pc, 
though these numbers are not strongly constrained. The 
X 2 = 1 surface spans R = 50 - 140 pc and a = 0.2 - 0.8. 
Moreover, the assumption of a virial parameter equal to 
1 (i.e., that Equation [T41 holds) is questionable both be- 
cause we neglect support by magnetic fields, non-circular 
geometries, and surface pressure terms (while considering 
substructure inside of a larger cloud) . Even more gener- 
ally, the fundamental assumption that clouds or parts of 
clouds are virialized is not certain to hold. 

Despite these concerns, Figure [8] does demonstrate 
that a simple model — CO emission nested inside a 
larger sphere of H2 — can relate dynamics measured from 
molecular line emission and H FIR . The best fit radius, 
R = 70 pc, is quite similar to that needed to achieve 
the extinction threshold for CO emission (Ay ~ 1) us- 
ing our adopted DGR and n k, 100 cm -3 — a typical 
average volume density for Galactic GMCs and perhaps 
appropriate for the diffuse gas between dense molecular 
clumps in the SMC. These three numbers combine to 
yield a depth of ~ 60 pc. Meanwhile, the density pro- 
file is similar to the a = 1 common ly used to describe 
Galactic clouds (|Solomon et al.lll987t ). 

The strong dependence of M v i r /M F 2 R on the size- 
scale sampled at least partially motivates the discrep- 
ancy between CO-to-H2 conversion factors measured us- 
ing CO observations and those derived from dust. At 
the high resolutions achieved by millimeter-wave inter- 
ferometers in Local Group galaxies, CO-emitting clouds 
are resolved from their surroundings. By concentrat- 
ing on these clouds, one samples only dense regions 
where CO is well-shielded by dust. This naturally leads 
to relatively modest conversion factors. On the other 
hand, dust measurements and dynamical measurements 
made on larger scales sample the whole complex. In 
the SMC this appears to includes a large amount of 
poorly-shielded gas and such methods therefore return 
significantly larger conversion factors. One manifestation 
of this phenomenon is that dynamical mass determina- 
tions from CO measurements with larger physical beam 
sizes often return systematically and significantly higher 
conversion factors than those obtai ned from CO mea- 
surements in much smaller beams (iRubio et al.l fl993b; 
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Fig. 8. — Reconciling dynamics and H!p R in N83. (top left) The ratio of virial mass, M v - iv , to total H2 mass, Af^j i (j/-axis) expected 
for the simple case where M v i r traces only an inner portion of a cloud (the fraction traced is shown on the rr-axis). Each line shows a cloud 
with a different density profile, (top right) N83 divided into concentric regions denned by CO intensity. We measure M v i r and M^i^ for 
each region, (bottom left) Results of fitting the models in the top left panel to M^ Y /M^^ measured from the regions in the top right 
panel. The x-axis show the power law index of the cloud density profile; the j/-axis shows the cloud radius. Contours show reduced \ 2 , 
starting at 0.5 and increasing by a factor of 2 each step. The white cross marks the best-fit model (p oc r -0,6 , R = 70 pc). (bottom right) 
M v i r /M^2 R as a function of region radius (black points), along with the best fit model (gray line). 

IWilsonl[l995l : lIsTael l2000t IBolatto et aT] 12001) . For in- 
terferometer measurements to properly sample the full 



cloud structure a multi-scale analysis, such as that pre- 
sented here or the mo re rigorous "dendogram" approach 
recently described by iRosolowskv et alj (|2008l ) , is neces- 
sary. 

Although our dynamical and dust-based results appear 
consistent with this simple picture, other recent results 
suggest a more complex relationship between the two 
measurements. Bot et al. (2009, in prep.) recently 
measured the relationship between sub-millimeter dust 
emission and CO-based dynamical masses in the south- 
west part of the SMC Bar. Even after controlling for 
contamination by an extended superstructure of CO-free 
H2, they find that virial masses are systematically lower 
than dust-based H2 masses on the scale of individual 



CO-bright regions. This might arise if clouds are short- 
lived (i.e., presently collapsing) or partially supported by 
magnetic fields. Alternatively it may reflect altered dust 
properties in dense cloud cores. The virial-dust discrep- 
ancy measured by Bot et al. and the multiscale virial- 
dust measurements presented here can both be readily 
applied to simulated clouds and multi-tracer observations 
of Galactic GMCs. It will be interesting to see whether 
these measurements can be replicated purely by altering 
the CO-emitting surface inside of a cloud (as it appears 
from our simple model) or if they constrain SMC cloud 
structure to be genuinely different from that in the Milky 
Way (as appears to be the case from the Bot et al. re- 
sults). 

6. SUMMARY AND DISCUSSION 



Structure Of A Low Metallicity CMC 



13 



We combine far infrared emission, CO line emission, 
and a 21-cm H I map to study the structure of CO, dust, 
and H2 in the SMC star forming complex N83. 

Two recent surve ys of the SMC using Spitzer (S 3 MC 
iBolatto et aflUM and SAGE-SMC, Gordon et al., in 
prep.) allow us to estimate the distributions of dust and 
H 2 at high spatial resolution. We calibrate a method to 
derive the equilibrium dust temperature, Td us t, and op- 
tical depth at 160/im, r 160 , along the line of sight using 
only Spitzer data. Applying this method and assum- 
ing that the diffuse ISM of the SMC Wing is mostly 
H I, we determine the dust-to-gas ratio (DGR) using the 
T160 and H I maps. We find T160 to be a good tracer of 
A(H I) with Ti 60 = 1.4±g;| x 1(T 26 cm 2 N (H), implying 
a DGR 17^g° (la) times lower than that in the Solar 
Neighborhood. High residuals about the T160 - iV(H I) 
relation come almost exclusively from regions of active 
star formation, with the largest residuals from N83 it- 
self. The most likely origin for these high residuals is 
dust associated with H2, though several important sys- 
tematic uncertainties remain unconstrained (Appendix 
lAl . Considering several pieces of evidence (the metal- 
licity of the N84C H II region, UV spectra of a nearby 
star, and the DGR in nearby diffuse ISM) we adopt 
a DGR of N(R)/E(B - V) « 5 x 10 22 cm" 2 mag" 1 
(tj.60 = 2.8 x 1(T 26 cm 2 A(H)) for N83 itself, but note 
this as a significant uncertainty with the plausible range 
spanning N(H)/E(B - V) = 3-10 x 10 22 cm" 2 mag" 1 . 
Combining this DGR with t 160 and the measured H I 
distribution, we derive a map of Hf 11 * in N83. 

Comparing CO intensity, kinematics, dust, and H2 we 
find: 

1. The CO-to-H2 conversion factor averaged over the 
part of the N83/N84 region mapped by SEST 
is very high, 4-11 x 10 21 cm -2 (K km s -1 ) -1 or 
« 20-55 times the Galactic value. Despite the 
large discrepancy from the Galactic A co , there is 
reasonable agreement between the distributions of 
CO and H2 traced by dust: a rank correlation co- 
efficient ~ 0.7 relates the two over the SEST field. 

2. Bright CO is more confined than H2, so that Aco 
varies across the region, with the lowest (most 
nearly Galactic) values near the CO peaks. The 
magnitude (or existence) of an extended, truly CO- 
free envelope is a sensitive function of the adopted 
DGR. Our best estimate is that such an enve- 
lope does exist, with Eh2 ~ 100 M Q pc~ 2 where 
Ico ~ 0. 

3. CO emission is a function of line-of-sight extinc- 
tion, which we estimate from t 160 . Bright CO emis- 
sion is largely confined to regions with Ay > 2 mag. 
This agrees we l l with modeling of SMC clouds by 
iLequeux et al.l (|1994l ) and roughly matches what 
is seen in the Milky Way. This result is robust to 
most of the systematic uncertainties that affect our 
determination of H2 . 

4. A simple model can reconcile dynamical masses 
(measured from CO) with H2 (measured from 
dust). In this model, CO emission comes a surface 
within the cloud while dust emission traces all H2 
along the line of sight. The best-fit density profile 



and radius are p oc r~ o e and R = 70 pc. These are 
not strongly constrained, but the density profile is 
similar to that inferred for Galactic clouds and the 
radius is consistent with that required to achieve 
Ay « 1 mag for our adopted DGR and a typical 
molecular cloud density. 

These results — particularly the confinement of in- 
tense CO to regions of relatively high line-of-sight ex- 
tinction — are all consistent with the selective pho- 
todis sociation of CO relative to H2 at low met all icities 
(e.g-lMalonev h Blacldll988t iRubio et al.lll993bl lal: llsraell 
1997b; lBolatto et al.lll999t ) In this scenario, the distri- 
bution of CO emission is largely driven by need for dust 
to shield CO from dissociating radiation. The underlying 
distribution of H2, while subject to significant systematic 
uncertainties, appears similar to that in a Galactic GMC 
complex. 

If the distribution of CO emission is indeed largely de- 
termined by dust shielding, then we expect that the ra- 
tio of CO emission to H2 mass will depend sensitively 
on both the local DGR and the radiation field incident 
on the cloud. These effects may largely cancel in more 
massive spiral galaxies, yielding a CO-to-H2 conver sion 
factor that is fairly robust (e.g., Wo lfire et al.lll993[ ). In 
low-mass galaxies, which have high radiation fields and 
low DGR, they will tend to compound, producing ex- 
tended envelopes of H 2 with little or no associated CO. 



From recent large surveys of the Magellanic 
Clouds at i nfrare d and millimeter wa v elengths (e.g.. 



Fukui et al.l 119991: 


Mizuno et al.l 120011: iMeixner et all 


2006: IBolatto et al. 


120071: lOtt et alJ 120081 Gordon et 


al., in prep.), it will be possible in the next few years 



to fill the right panel in Figure [7] with points from 
across the Clouds. This will allow the quantification of 
the radiation field (and perhaps density) as a "second 
parameter" in the Iqq-Av relation. It may also allow 
an improved calibration of Aco as a function of both 
DGR a nd loca l radiat ion field, extending the pioneering 
work bv llsraell (|1997bl ) to the scale of individual clouds. 

Even with such data, it is unclear if CO emission can 
remain an effective tracer of H2 on the scale of individual 
clouds. Tracing local variations in DGR and radiation 
field to apply a spatially variable Aco may not be pos- 
sible or practical. Of course, CO is already well-known 
to be a flawed tra cer of H 2 within Galactic clouds (e.g., 
iPineda et al.ll2008l ) but retains significant utility for trac- 
ing H 2 on large scales. Over a sizable portion of a galaxy, 
variations in the radiation field and DGR may average 
out and allow a calibration to work at a basic level. Given 
that the options to trace H2 in low-metallicity galaxies 
remain limited, a combination of dust and molecular line 
emission is likely to be the only widely available option 
in the near future. Herschel spectroscopy of the [CII] line 
and Fermi observations of 7 ray emission from the Mag- 
ellanic Clouds, while both likely to illuminate the issue 
significantly, will only target a small sample of galaxies. 

We thank the anonymous referee for a detailed and 
helpful critique. We thank Henrik Beuther for helpful 
comments on a draft of the paper and Mark Krumholz for 
a helpful discussion. We acknowledge the use of NASA's 
Astrophysics Data System Bibliographic Services. 
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Fig. 9. — (left) The average H I spectrum over the region of high residuals (black) and spectra from individual lines of sight in this area 
(gray). The spectrum of CO emission (with an arbitrary normalization) is shown below the H I. (right) The distribution of opacities in the 
(integrated) 21cm line required to explain the residuals in highest contour in Figure [4] Although individual spectra show some evidence 
of optical thickness, we see no clear signature of self absorption. The line-integ rated values of T2i C m required to explain ri60 in N83 are 
mostly higher than the peak values of T2i cm measured anywhere in the SMC by[Dickcv ct al. (200(3). 
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APPENDIX 

SYSTEMATIC UNCERTAINTIES IN N (Hf IH ) 

Several systematic uncertainties may affect N (Hf IR ) but are hard to quantify and so not reflected in our Monte 
Carlo estimate of the uncertainties. Here we discuss these for the specific case of N83 (for a more general discussion 



16 



Leroy et al. 



□ 
S 



3 
I 



-TB'ED'D 1 




B 100 



□i b iimr 



13=SQ" 




a 4 

In, B— L [K km B" 1 ] 



Fig. 10. — (left) Ho emission (gray scale) with S^ 1 ^ shown in contour (both at 36" resolution). Although Ha and Hj 111 roughly coincide 
on large scales (Figure 0), the detailed distributions are not a good match, (right) The effect of changing DGR to the most extreme 
plausible values on the relationship between Sg2 R (y-axis) and Ico (x-axis, only J = 2 — » 1 shown). The black points show the lowest 
plausible DGR in N83, that found in the nearby diffuse gas. The gray points show the highest plausible DGR, ~ 3 times this value. 



see Ilsraelll997al) . We find no strong reason to doubt that Equation [7] yields an approximate estimate of iV(H2). N83 
appears unlikely to harbor a significant population of cold dust and we do not observe compelling evidence that dust 
traces mostly warm ionized gas or high optical depth H I. There is likely some blending of populations along the 
line of sight, but the magnitude of the effect is unclear. Grain processing is largely unconstrained, but we note the 
dissimilarity between N83 and the dense, cold cores where these effects are usually discussed. 

Blending of Populations Along the Line of Sight: N83 is a dense, active region and the line-of-sight distance through 
the SMC may be very long. As a result, the observed dust emission may represent a blend of several dust populations 
with different temperatures. The likely effect is that we overestimate the a verage Tdust along the line of sight and 
therefore underestimate ti 6 o and H™ (e.g., see tests on simulated clouds bv lSchnee et al.|[2006[ ). 

Cold Dust: A related concern is that our longest wavelength data are at 160/im. As a result, we would miss any 
population of cold dust. In the Milky Way, when cold, molecular filaments can be isolated from embedded star 
form ation, they are often observed to have low dust temperature s (T < 15 K) and little out-of-equilibrium emission 
(e.g.. lLaureiis et aflll991l : lBernard et al.l[T999HStepnik et al.ll2003l ). As with blending of several populations, cold dust 
is most likely to be associated with the dense, molecular environment of N83. Missing cold dust would lead us to 
underestimate t\qq and A r (Hf IR ). 

Given the high Tdust in N83 and the presence of ongoing, vigorous star formation we consider it unlikely that there 
is a significant a mount of cold dust present. We attempt a simple test that reveals the presence of cold filaments in 
Galactic GMCs (jAbergel et al.l Il99^ iBoulanger et aT1ll998[ k we take the median Iiea/ho over the region, scale the 
70/im map by this value, and subtract it from the 160/im map. This should reveal the location of any local 160/im 
excess, a likely signature of cold dust. We find no such excess associated with N83 as a whole or the CO peaks in 
particular. 

Other Gas Phases: We refer to the results of Equation [7] as " IR " but this is actually an estimate of all gas not 
traced by the 21-cm transition. Some of this might be high optical depth H I or warm ionized gas. Neither appears to 
be a plausible explanation for the majority of such gas in N83. This agrees with results from the Milky W ay, where 
excess dust emission identified in a similar way also appears to correspond mostly to H9 (|Dame etmi200l . 

The right panel in Figure [9] shows the H I opacities required to acc ount for T\m i n N83 given our adopted DGR. 
These values, T2i cm — 2-4, are higher than those implied by the fit of iDickev et al.l (|2000h . which yields a maximum 
line-integrated T2i cm ~ 0.55 (correction factor ~ 1.3) near N83. Indeed, most of the line-integrated values of T2i cm in 
Figure El are higher than any of the peak T2i cm values (i.e., T2i cm in the most opaque velocity channel) measured by 
IDickev et al. (2000) in the SMC (maximum ~ 1.7). though that s tudy did not probe any star-forming peaks; toward 
the starburst region 30 Doradus in the LMC IDickev et al.l (|1994j ) found peak T2i C m values of ~ 2, which is still too 
small to achieve the line- integrated value of T2i cm required account for t\qq in N83. The 21cm spectra do show some 
evidence of optical thickness at a brightness temperature of ~ 120 K, but no clear signs of self-absorption at the 
velocity of the CO peak (left panel in Figured]). We cannot rule out unlucky geometry, but achieving line-integrated 
optical depths of 2-4 without invoking a contrived scenario appears difficult. 

Warm ionized gas also seems unlikely to account for most of H!p R . The left panel in Figure [TO] shows contours of 
EES 1 over an Ha image (at matched resolution) in the SEST field. Although high r 160 residuals correspond to Ha 
emission on large scales, the detailed distribution is not a particularly good match. The rank correlation coefficient 
between Ha and Hf IR over the area observed by SEST is ~ 0.1, much lower than the 0.7 relating Hf IR and CO. Ha 
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emission is proportional to J n 2 dl and so obviously a flawed tracer of the true warm gas column (J ndl), but the poor 
correspondence on small scales still argues that most Hj 111 is not actually warm ionized gas. 

Dust Processing in Molecular Clouds: A significant but hard-to-constrain uncertainty in Equation[7]is if and how dust 
properties vary between N83 and the surrounding ISM. The most likely variations are increases in the FIR emissivity 
or the DGR. In the Milky Way, the FIR emissivit y of dust (tfirMv) does appear to increas e towards dense regions , 
increasing by ~ 30-50 % above Ay ~ 1 mag (e.g., lArce fc Goodman] 19991 : iDutra et al.ll2003t ICambresv et al.l 120051) . 
ICambresv et al.l (|2001(l . IStepnik et al.l (f2003h . a nd ICambresv et al.l (|2005l ) argue that this is due to the creation of 
fluffy dust grains with low albedos ( Dweklll997f) via grain-grain coagulation or accretion of gas. At the same time, 
build -up of existing grains in molecular clouds and dust creation in Type II supernovae or stellar winds (e.g.. [Dwekl 
1998) may cause the DGR ratio near star-forming regions to be higher than in the surrounding ISM. 

The magnitude of grain growth in GMCs remains very poorly constrained and in an active environment like N83 
it will be balanced against grain destruction (e.g., in shocks). Further, the high dust temperatures, low integrated 
extinctions {Ay < 2 mag almost everywhere), and weak CO emission in N83 are a far cry from the high density, high 
extin ction environments in which grain coagulation or the form ation of icy ma n tles a re usually modeled or observed 
(e.g.. lOsscnko pTfc Hennin"glll994f ). Moreover, as pointed out bv lBernard et all (l2008f). increase d emissivity in Milky 
Way clouds is often associated with diminished small grain emission (e.g.. ISchnee et all 12008( 1. while N83 exhibits 
increased Im/Iim compared to its surroundings. 

Because it is unclear what, if any, grain processing is at work in N83, we make no correction to the emissivity. If 
dust in N83 indeed has a high emissivity compared to the diffuse ISM, we will derive values of both iV(H2 IR ) and Ay 
that are too high. Note that our adopted DGR is already twice that in the surrounding diffuse gas. Increasing or 
decreasing the adopted DGR will not affect Ti6o or Ay, but will lower or raise iV(Hf IR ). 

Effect of Changing DGR on the CO-H2 relation: The exact value of the DGR in N83 is the largest systematic 
uncertainty in our analysis. We discuss the constraints on this quantity in £ 14.31 In the right panel of Figure HQl we 
illustrate the relationship between H2 and DGR in the limiting cases: DGR equal to that in the diffuse ISM of the 
SMC Wing (black) and DGR equal to three times this value (gray). 

There are two main conclusions to draw from this comparison. First, the existence and magnitude of a truly CO- 
free H 2 envelope (the y-intercept of the points) depends sensitively on the adopted DGR; the lowest plausible value 
is (partially by construction) consistent with no envelope and the highest value implies an envelope with surface 
densities ~ 200-400 M Q pc -2 , 1-2 times the average surface density of a Galactic GMC. Second, the average CO-to- 
H2 conversion factor varies between 10 and 100 times Galactic over the full range of possible DGR. The qualitative 
behavior (meaning the presence of bright Iqo only above a certain Ejj2 threshold) remains the same. We emphasize 
that the relationship between Ay (or ri6o) and Ico is unaffected by the choice of DGR. 



